Critical hybridization for the Kondo resonance in gapless systems 
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We study the Kondo resonance in a spin- 1/2 single impurity Anderson model with a gapless 
conduction band using the equation of motion approach in order to obtain the impurity spectral 
function. We study two different scenarios for gapless systems: a purely power-law energy depen- 
dence of the density of states or a constant density of states with a gapless behavior near the Fermi 
level. We demonstrate that strong electron-electron correlations lead to a sharp peak in the impu- 
rity spectral function in the case of a large exchange coupling (J > J c ) or equivalently, a strong 
hybridization (V > V c ). This Kondo-like peak emerges much below the Fermi level in the case of a 
strongly depleted density of states. These results are compared with the ones from renormalization 
group approaches. 

PACS numbers: 



I. INTRODUCTION 



The advent of graphene has renewed the interest in 
the Kondo effect in gapless systems. [l| It is known that a 
strong depletion of the density of states of a conduction 
band near the Fermi level modifies the usual behavior 
of magnetic and nonmagnetic impurities. @, i Impuri- 
ties in gapless systems have been addressed, for instance, 
in the context of d-wave superconductors 0, [f| @ and 
more recently graphene, Q using a variety of methods, 
from renormalization group approaches [g, d, OH Hl| to 
the large N analysis [12| and T-matrix calculations, [l 3| 
The most relevant modification in gapless systems of 
the usual Kondo behavior is the requirement of a fi- 
nite critical Kondo coupling strength in order to observe 
the Kondo resonance. [2|, |7| While several renormaliza- 
tion group studies have addressed this point [1, H, 0] , an 
equivalent study from the point of view of the impurity 
spectral function in the case of the single impurity An- 
derson model has not been carried out as far as we know. 
In the present manuscript, we address this problem using 
the equation of motion (EOM) approach which has been 
successful in describing the Kondo resonance in systems 
with a finite density of states at the Fermi level. [lj] The 
nonmagnetic case when the on-site Coulomb repulsion is 
zero (U = 0), and the strong on-site Coulomb repulsion 
limit (U — * oo) of the gapless single impurity Anderson 
model are considered. Note that above the Kondo tem- 
perature, the U — > oo spectral function should recover 
the U = features. We will show that even for U = 0, 
an additional peak may appear in the spectral function 
due to the depletion of the density of states which should 
not be confused with a Kondo peak. We will demonstrate 
that the critical exchange coupling depends strongly on 
the density of states profile in the gapless region. Further- 
more, the strong electron-electron correlations can pro- 
duce a strong Kondo-like peak in the impurity spectral 
function even much below the Fermi level in the depletion 
region. This Kondo-like peak emerges if the hybridiza- 
tion between the conduction band and the impurity is 



larger than a critical value, which in turn is determined 
by the energy dependence of the density of states in the 
depletion region, and grows with decreasing temperature. 



II. GAPLESS ANDERSON HAMILTONIAN 

The single impurity Anderson model is given by 
H = e k n krT + 2J E n dcr 

ka a 

-^U^h da h d5 + V^{cl d a + d+c ka ). (1) 



ka 



where the standard notations[l5| are used: V is the hy- 
bridization between conduction and impurity states, and 
E is the impurity energy level. In order to describe gap- 
less systems, a density of states with a power-law energy 
dependence is usually considered in a certain energy in- 
terval around the Fermi energy which we set at £f = 0. 
We will assume that the density of states of the conduc- 
tion band is constant outside this interval, p{e) = poj for 
D 1 < \s\ < D, and p(e) = pi(e) = a|e| 7 for |e| < D x with 
a = po/Dj. Therefore, p(e) is zero at the Fermi level. 
The bandwidth is 2D. In Fig. [TJ the density of states 
for several values of 7 is shown. In graphene, 7 = 1. If 
D\ = 0, then one has a constant density of states with 
the finite bandwidth. If D\ = D, only the power-law 
dependence of the density of states is considered. 7 = 00 
corresponds to a system with a 2Di gap. Below, we will 
discuss in detail the case Eq < —D\. 



III. EQUATION OF MOTION METHOD 

The equation of motion approach is basically the suc- 
cessive application of the result wGab(w) = ({A, B}) + 
(([A,H];B)) U to the impurity Green's function and to 
the other Green's function generated in the process. [f3| 
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FIG. 1: The density of states p(uu) of the conduction band 
used in this article. The density of states is constant in the 
energy range Di < \lu\ < D and has a power- law dependence 
p(lu) — a|tj| 7 for \lu\ < Di. The density of states profile for 
several 7 values is displayed. Note that 7 = 00 implies the 
existence of a 2D± energy gap. 



Here we have adopted the Zubarev notation for the re- 
tarded Green's function, Ga.b{u) = ((A; -B)) w .[l6| As 
usual, it is implicit that lu — > u> + irj where 77 is a in- 
finitesimal positive constant. In the case of the Anderson 
model, one obtains the following EOM for the impurity 
Green's function 



v- 



tO - Ek 



((d a ;dt)) u = l 

+ U{(h d ad a ;d+)) U . (2) 



If the on-site repulsion between electrons is neglected 
(U = 0), this equation becomes a closed EOM. Below, 
we discuss this {7 = case in detail before addressing 
the U — > 00 limit for two reasons. First, due to the de- 
pletion of the density of states at the Fermi level, the 
U = impurity spectral function will have a non trivial 
dependence on the impurity energy level and on the hy- 
bridization value V. In particular, it may display more 
structure than the typical Lorentzian profile. This more 
complex profile is the background against which the cor- 
relation effects will be observed in the U — * 00 limit. 
Second, the U — features should be recovered in the 
U — > 00 spectral function above the Kondo temperature. 



IV. RESONANT LEVEL (U = 0) 

The self-energy for U — can be easily obtained for 
integer 7 and is given by 



V 2 



LU - E k 



A(w) - iA(w), 



(3) 



FIG. 2: The real part \(u) of the self-energy E , Eq. |J5J, for 
the resonant level ((7 = 0) exhibits logarithmic divergences 
at the band edges and a non-monotonous behavior in the 
gapless frequency region. The dashed line represents the line 
(u) — Eq)/V 2 . The black points are the roots of lu — Eo — A(u;). 
For 7 = 00 (the gapped spectrum), additional logarithmic di- 
vergences are present at the gapless region edges. Parameters: 
Eo = -5, po = 1, V = 4.4, D 1 = 4 and D = 40. 



with 



AM 
V 2 

AM 
V 2 



t(lu) = ttpM + 77, 

am = bm + p„m 

{D + uu){D l - lu) 



po In 



(D-u){Di+u) 



(4) 



(5) 



where 



B(lu) = < 



au) 1 In 



ulu 1 In 



LU 



D 2 -u 2 

Dt+LU 



- lu 



7 odd, 



7 even, 



(6) 



and P n (u>) is a polynomial of degree n < 7. For 7 = 1 
and 7 = 00 (a gapped conduction band), P n (uS) = 0. In 
Fig. [21 the real part of the self-energy So is displayed for 
7 = 0, 7 = 1 and 7 = 00 (the gapped spectrum). 

The spectral function A(lu) = — 2 Im ((d a \ d^)) u for a 
system with finite bandwidth (gapless or not) is given by 



A(lu) 



2A(w) 



[lu - E 
2n6(D - 



A(w)] 2 + A 2 M 

D^li-A'Mi- 1 ^ 



(7) 



where lui are the roots of lu — Eq — A(lu) = 0. Note that 
the Dirac delta points outside the band continuum play 
an important role if the hybridization is not small. These 
points are essential in order to satisfy the sum rule of the 
spectral function. This can be more clearly understood 
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FIG. 3: The spectral function A(u), from Eq. (0, for the 
resonant level ((7 = 0) for 7 = 1 and several values of the 
impurity level energy: a) Eq = —15; b) Eo = —10; c) Eo = 
— 5; d) Eo = 0. Inset: The real part of the self-energy A(w) 
and the line (to — Eo)/V 2 . It is the non-monotonous behavior 
of A(w) that generates additional structure in the spectral 
function for certain intervals of the parameters. For small V , 



this structure disappears. Other parameters: p 
D 1 = 4 and D = 40. 



if the continuous term is written as 
1 2j(oj) 



1/2 



V 



1.4, 



V 2 



V 2 



El 
V 2 



AH 



(8) 



Therefore, for large V and a constant density of states 
most of the spectral weight is associated with the Dirac 
delta points outside the band continuum since the total 
spectral weight for |w| < D decays as 1/V 2 . As can be 
inferred from Fig. [2J for large V [or equivalently a small 
slope of the straight line (u> — Eo)/V 2 ], the Dirac delta 
points are located at large u>. For large w, A(w) ~ V 2 C /u> 
where C ~ pqD and therefore a;, = ±VVC leading to 
|1 — A'(cj)| _1 ~ 1/2. Thus the sum rule is satisfied. For 
small V, these Dirac delta points are located near the 
boundaries of the band due to the log-like behavior of 
A(u>). The same conclusion holds for a power-law density 
of states. 

For a constant density of states, X(uj) is monotonous 
for \u>\ < D and the minima of the first term in the 
denominator of the spectral function given by Eq. ([5]) 
are typically roots of u> — Eq — A(u>). In the case of a 
power-law density of states, additional minima appear 
due to the non-monotonous behavior of A(ui) within the 
gapless region. These minima are not necessarily roots 
of oj — Eq — A(u>), but lead to an additional structure in 
the spectral function. In Fig. [3l the spectral function for 
U = at several values of Eq is displayed. For general 7, 
one has a |w| 7 behavior of the spectral function near the 
Fermi energy, except for Eq = as can be observed in 
Fig. [3j It is worthwhile to emphasize that the impurity 
spectral function exhibits the same depletion at the Fermi 
level as the density of states of the conduction band. 

Note that the value of V determines the slope of the 



line (uj — Eq)/V 2 leading to different behavior of the spec- 
tral function. A small V implies a large slope and typ- 
ically only one root of u> — Eq — A(u>) for \co\ < D. For 
large V, one has an almost horizontal line intersecting 
the A(w) curve in Fig. [2j The same reasoning applies in 
the interacting case that will be considered below. Since 
the V 2 terms appear always associated to a pq factor in 
the spectral function, these parameters play equivalent 
roles. 

It is important to emphasize that a root of [u> — E — 
A(w)] is not a maximum of the spectral function if this 
root is within the gapless region of the density of states. 
In fact, given a function g{x) — x/[f 2 (x) + x 2 \ where 
f(x) is a function with a local peak at x = xq, then 
the corresponding peak in g(x) will be at x' = Xq{1 + 
1/[/'(il-o)] 2 }~ 1/2 and therefore x w x if /' > 1 or iwO 
if /' <C 1. Furthermore, if the root is near the Fermi level, 
the respective peak in the spectral function becomes very 
narrow, but the total spectral weight associated with the 
peak becomes very low. This can be easily concluded 
since the spectral function [Eq. ©] for 7 = 1 can be 
written in the neighborhood of the roots (if \u>i\ < D±) 
as c|oj| J-"(oj), where T{u>) is a normalized Lorentzian with 
center at uii/[l + K 2 ] and half-width \K/(l + K 2 )-uji\, where 
k = naV 2 /\l - A'(wi)| and c = n\l - A'(cj ? )|- 1 |k/(1 + 
k 2 )\. For small \u>i\, A'(uii) becomes large and k becomes 
small, and therefore one has a narrow peak with a low 
spectral weight. The linear factor leads to a shift of the 
maximum of the Lorentzian away from zero frequency to 

^/[l + K 2 ] 1 / 2 . 



V. THE KONDO RESONANCE (U = 00) 

In the U case, the equation of motion for the 
last term in Eq. ^ has to be considered. Adopting the 
Appclbaum and Penn approximation, [13] one finds the 
following expression for the impurity Green's function 
for U = 00, 



(w - Eq - So - Ei - E 2 ) ((da-, = 



LO - Sk 



where 



E = y2y (4^fcq) 
W — £k 

k,k> 



= ^ 3 E 



(d£c ks ) 



k.k 



f (lo - e k ){u) - EfeO ' 



(9) 

(10) 
(11) 



and E is given by Eq. (J3J). In order to describe qual- 
itatively the Kondo resonance, one has to understand 
the temperature behavior of the low order self-energy 
Ei. For sufficiently small V, one can neglect E 2 and 
the Kondo resonance will be generated by Si, which is 
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FIG. 4: The real part of the zero temperature self-energy 
\i(uj) = Re [£i(cj, 0)]/V 2 for several values of 7. In the case 
of 7 = 00, there is a gap, N(e) = at \ui\ < Di, which leads 
to a logarithmic divergence of the real part of the self-energy 
at the gap edge. Other parameters: po = 1, D\ = 4, and 
D = 40 or D = Di. 



approximately given by 
Ei(w,T) „ 



mp(u})f{u}) 
dw'p(Lj')f(w')P 



(12) 



The Fermi-Dirac distribution function f(x) sets a high 
energy cutoff in the previous integral which, for 7 = 1 , is 
approximately given by 



Re 



Ei(w,T) 



V 2 



Po In 



D 



uio In 



- a{D 1 - k B T) 

,Jw 2 + (ak B T) 2 
\Dx+w\ 



(13) 



for k B T < Di and a ~ tt. We would like to note that in 
our numerical solution of Eq. ([9]), the self-energy E2 is 
taken into account. 

The previous expression should be compared with the 
Lacroix result for a constant density of states. 



Re 



1/2 



~Po In 



^ 2 + {nk B Tf 
D 



(14) 



The most relevant difference between the previous two 
expressions is the temperature behavior. For a constant 
density of states, Re(Ei) diverges at zero temperature as 
u> — > 0. For the gapless density of states, as the temper- 
ature goes to zero, Re(Ei) remains finite at the Fermi 
energy (as shown in Fig. |4]) with a peak in the energy 
range —D\ < to < 0. Note that for a strong depiction 
of the density of states this peak grows and approaches 
the gapless region edge, converging therefore to the log- 
arithmic divergence of the real part of the self-energy of 
an equivalent system with a 2D\ gap. In Fig. [4] the self- 
energy for a linear and the quartic density of states with 



FIG. 5: The spectral function A{u>), for U = 00, 7 = 1, and 
several values of the hybridization V. Curve a) displays the 
typical V < V c spectral function characterized by the ab- 
sence of a Kondo resonance. The Kondo peak appears as 

V is increased above a critical value V c [curve b)] and dis- 
appears above the Kondo temperature T c [curve c)]. Other 
parameters: Eo — —3.2, po = 10, D\ = 1, D = 20 and 
a) V = 0.12, k B T = 0.05; b) V = 0.24, k B T = 0.05; c) 

V = 0.24, k B T = 10. 



D = D\ is also displayed. The self-energy of the latter is 
strongly reduced in the —D\ < lu < region due to the 
the absence of the log-like contribution which resulted 
from the — D < u < —D\ constant term in the density 
of states. 



A. The critical hybridization 

The critical temperature for the Kondo resonance 
in the equation of motion approach [l4| is associated 
with a temperature below which a new solution of 
Re[((d CT ; = appears. This solution is due to 

an additional intersection between the line u — Eq and 
the real part of the self-energy Rc(Ei). In the case of 
a constant density of states, this solution occurs near 
lu = 0. Since the maximum value of the self-energy for 
the constant density of states (7 = 0) is 



Re 



£i(0, T) 



V 2 



-po In 



Tk B T 
D 



(15) 



and decreases with increasing temperature, above a criti- 
cal temperature one loses the additional intersection and 
the Kondo peak disappears. The condition E /V 2 ~ 
pahi(k B T c / D) determines the Kondo temperature. Note 
that this condition has a solution for any value of V, how- 
ever small it may be, due to the logarithmic dependence 
on temperature. 

In the case of the gapless density of states (with 7 > 1), 
the same reasoning can be followed in order to obtain 
the critical temperature. However, the real part of the 
sclf-cncrgy Ej docs not diverge at zero temperature for 
any frequency (excluding the band edges). Therefore, for 
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FIG. 6: The spectral function A(ui) for U = oo, 7 = 8, and 
several values of the impurity energy, En = —5, — 7 and —9. 
One can see that when Eq increases and approaches the edge 
of the depletion region (ui — —3, on this figure), the Kondo 
peak grows within the depletion region much below the Fermi 
energy. Other parameters: pn = 10, D\ = 3, D = 20, V = 
0.2, and k B T = 0.05. 



a sufficiently small hybridization, the Kondo peak does 
not occur. In the case of p(uj) = a|cj| 7 0(Z?i — \u>\), the 
real part of the self-energy for large 7 can be obtained 
expanding the density of states at u> = —D\ and is given 
approximately by 



Re 



Si(faJ,0) 
V 2 



aDj + aDj 



x In 



1 



1-7(1 + ^ 



1/7 



1 



(16) 



This function has a maximum at Lo max given by uj ma x + 
Di ~ 3£>i/47 and approximately of the value of p . 
The constant behavior of the density of states in the en- 
ergy interval ~D < oj < ~D\ leads to an additional 
P0V 2 In \ {D + lo)/(Di + oj)\ term in the real part of the 
self-energy, Re 0)]. For D 3> D±, the maximum of 

the real part of the self-energy is shifted to Lu max + D\ ~ 
D1/27 and its value is approximately poV 2 In (27D/D1) 
which is the same as the former maximum value except 
for an enhancement factor \n(2jD/Di). Note that this 
enhancement reflects the contribution of the constant 
density of states region. In fact, for large 7, the behav- 
ior becomes similar to a gapped system. In Fig. 2J the 
difference between the maximum values of Re[£i(w,0)] 
for D > D[ and D = D\ is clearly observed as well as 
the shift of the maximum position (with increasing 7) 
towards the negative edge of the gapless region. 

The critical hybridization is the minimum V which is 
necessary to have an intersection between the line u> — Eq 
and Re [Ex (a;, 0)] at zero temperature. For Eq < — D%, 
one has to take into account that the maximum occurs 
at ijJ max and the condition for the critical hybridization 
becomes for 7 > 1 



En 



V 2 



aDj \n{2 1 D/D 1 ) 
aDj, 



D = Dl 



(17) 



It is necessary to note that the density of states used 
in the renormalization group approach [2j corresponds 
to D = D\. For not so large 7, this density of states 
is reasonable, but in the 7 — ► 00 limit, this density of 
states becomes unrealistic since the electron states are 
concentrated in a narrow region at the band edges. Equa- 
tion (|17[) should be compared with the expression for 
the critical exchange coupling J c of the spin- 1/2 Kondo 
model with a gapless conduction band, J c ps j/aDj = 
7/Po-0 The Kondo model can be obtained from the An- 
derson model within perturbation theory fl~5j leading to 
an exchange coupling Jkk' ~ V 2 /(eu ~ Eq) which is an- 
tiferromagnetic for energies £k ~ 0. The resonant scat- 
tering of the conduction electrons in that energy region 
leads to the usual Kondo behavior with J ~ — V 2 /Eq. 
In fact, the Kondo limit corresponds to taking the limits 
V — > 00 and Eq — > —00 (as well as U — > 00) in the An- 
derson model. That is why the correction due to ui max is 
not observed in the renormalization studies of the Kondo 
model 0, 01 and has been also i gno red in renormalization 
studies of the Anderson model. |1C| 

For a gapped conduction band (7 = 00) with a 2Di 
gap, it is straightforward to conclude that the exchange 
coupling becomes J ~ V 2 j(D\ — Eq). In the case of a 
gapless conduction band, one also has to take into ac- 
count the depletion of states near the Fermi level and 
a correction in the previous expression should appear. 
In fact, a continuous evolution between the two cases 
should occur as 7 is varied. The real part of the self- 
energy, Eq. (|12|) . and in particular its maximum, reflect 
this evolution. Consequently, in the gapless case one has 
J ~ V 2 / {uj m ax — E Q ) and J^ 1 is given by Eq. (fT7)) which 
shows two different scenarios for the gapless energy re- 
gion in what concerns the Kondo resonance. For a pure 
power-law density of states with 7^1, the equivalent 
critical coupling will be J c ~ 1 /aDj = 1/po and there- 
fore is reduced in comparison with the renormalization 
group result 0] in the spin- 1/2 Kondo model by a fac- 
tor 7. For D 3> D\. the gapless energy region provides 
a cutoff to the log-like behavior of the real part of the 
self-energy for 7 = 00 and 

J" 1 - aDj ln(27Z?/£>i) = Po ln(2 7 £»/ J D 1 ) (18) 

and therefore an important logarithmic factor is present. 
Note that according to Eq. (jT5J) as 7 becomes larger, the 
critical coupling goes to zero, contrary to the renormal- 
ization group result, J c « j/po- We suggest that this dis- 
crepancy is due to the unrealistic behavior of the density 
of states used in the renormalization group approach, see 
the discussion above. The behavior described by Eq. (fTS)) 
reflects the existence of a Kondo-like peak in the impurity 
spectral function in the strong depletion limit, 7 3> 1 (see 
Fig. [5]). This phenomenon can be related to the Shiba 
intragap states that occur near the band edges in s-wave 
superconductors. [HJ 

In Fig. O the spectral function for U = 00 is displayed 
for several values of V and for temperatures above and 
below the Kondo temperature. These results were ob- 
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FIG. 7: Phase diagram of the gapless Anderson model in the 
7 — J^ 1 plane. The results presented in this paper are de- 
scribed by the solid red line. Above the solid red line (region 
I), i.e., in the region with a sufficiently small exchange cou- 
pling (or equivalently, a small hybridization parameter V), 
the Kondo-like peak is not formed. Strong electron-electron 
correlations lead to a strong peak in the impurity spectral 
function in region II, i.e., the region of a strong exchange 
coupling (J > J c ) and a strong hybridization (V > V c ). The 
blue dotted curve represents the renormalization group result 
of Withoff and Fradkin @]. 



tained numerically from Eq. ([9]) taking into account all 
self-energy terms. Curve a) displays the typical behavior 
of the impurity spectral function when (u) max — Eq)/V 2 
is above the critical value and no Kondo resonance is ob- 
served (even at zero temperature). In curve b) of Fig.O 
a Kondo peak is present since (uj max — Eq)/V 2 is below 
the critical value. Note that this Kondo peak does not 
occur at uj = reflecting the fact that the maximum of 
the real part of the self-energy Si occurs at uj max in con- 
trast to the usual Kondo behavior. A broader and lower 
Lorentzian profile is also observed in curve b) in Fig. [5] 
reflecting a larger value of the hybridization. Increasing 
temperature above the Kondo temperature, the narrow 
Kondo peak disappears but a small broad peak may re- 
main due to the depletion of the density of states at the 
Fermi level as explained in the discussion of the U = 
case. This is the situation displayed in curve c) of Fig. [5] 
For large 7, the Kondo peak emerges within the de- 
pletion region near its negative edge as shown in Fig. [5] 
Lowering the impurity energy level leads to the disap- 
pearance of the Kondo peak in agreement with the con- 
dition given by Eq. (fT7|) . Note that the Lorentzian pro- 
files have similar heights and widths in Fig. [5] since the 



hybridization energy is the same for all curves. 

The critical coupling relation [Eq. (fTBll ] leads to the 
qualitative phase diagram of Fig. [Jj The renormalization 
group result is also displayed in this figure. A consid- 
erably larger region (region II) of resonant Kondo-like 
correlations induced by the strong on-site Coulomb re- 
pulsion is present in comparison with the renormalization 
group result (dotted blue curve in Fig.[7J). This deviation 
becomes larger as 7 grows. Note that, in this paper, we 
have not addressed the 7 < 1 case. However, it is reason- 
able that as 7 goes to zero, the constant density of states 
behavior is recovered. Our belief is that the renormaliza- 
tion group approach describes correctly the 7^0 limit 
and therefore, J^ 1 cx 1/7 for 7 -C 1. The dashed part of 
the red curve in Fig. [7] reflects this assumption. 



VI. CONCLUSION 

In conclusion, we have determined the critical hy- 
bridization for the Kondo effect in the gapless Anderson 
model. Previous works have used renormalization group 
approaches to the Kondo effect in gapless systems and 
considered a purely power-law density of states. Here we 
have shown that if one considers a more realistic density 
of states, the critical coupling expression is reduced by a 
logarithmic factor that reflects the fact that the Kondo 
effect is dominated by the contributions of the density of 
states away from the gapless energy region. We showed 
that in a gapless system, strong electron-electron corre- 
lations between conduction and impurity electrons result 
in a sharp Kondo-like peak if the hybridization is larger 
than a critical value. The critical hybridization depends 
on the behavior of the density of states in the deple- 
tion region and on the energy of the impurity level. The 
low energy critical exchange coupling J c for the emer- 
gence of the Kondo resonance reflects the depletion of 
states at the Fermi energy and therefore is determined 
by the contribution of the states near the gapless region 
edge. We have demonstrated that in gapless systems, due 
to the strong on-site Coulomb repulsion, a sharp Kondo 
peak emerges much below the Fermi energy if the den- 
sity of states is strongly depleted. This Kondo peak has a 
strong dependence on temperature. Although the results 
presented above rely in a specific form of the density of 
states, we believe that any system with a well behaved 
density of states (except for its gapless behavior at the 
Fermi energy) should display similar behavior. 
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